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During recent years several very thorough and interesting 
treatises have been devoted to the investigation of systems of 
spars and ribs. These researches represent considerable prog- 
ress, especially for the statics of cantilever wings. From 
them we have obtained a clearer idea of the effect of the ribs 
and have also been enabled, in a few special cases, to calcu- 
late, by relatively simple methods, the magnitude of the prin- 
cipal effect of the unions of the spars and ribs, namely, the 
reduction of the bending moment at the fixed ends of the spars. 

The present paper treats less of the effect of the union 
than of its influence on the torsional rigidity of the wing. 
The calculations are carried out for a two-spar wing of con- 
stant cross section, in which the ribs are replaced by a con- 
tinuous member of constant rigidity (assumption of an infinite 
number of ribs). The bending rigidity of the ribs is assumed 
to be very great in comparison with their torsional rigidity, 
as is the case in actual practice. 

♦"Ueber die Torsionssteif igkeit eines freitragenden Flugels mit 
konstantem Holm- und Rippenquerschnitt " from Luftf ahrtforschung, 
June 26, 1928, pp. 79-90. 
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The case under consideration is strongly idealized. In 
reality the bending and torsional rigidities of the spars vary 
throughout the length of the wing* Nevertheless, the accurate 
numerical evaluation of this idealized case is of some inter- 
est, because it gives us an accurate idea of how the ratio of 
the bending to the torsional rigidity of the spars affects 
the torsional rigidity of the wing as a whole. Above all, it 
enables us to tell when the torsion of the whole wing may be 
regarded merely as the result of the various flexures of the 
spars and when, on the other hand, the wing may be regarded 
simply as a twisted girder. In this sense the designer can 
make approximate calculations, even in the most complicated 
cases, on the basis of the following method. 

Before proceeding to the calculation, we will take a look 
at the reports which have already appeared, especially during 
the preparation of the present paper. 

1. We find an article by Ballenstedt . * He discusses a 
wing structure with two like spars of constant cross section, 
and calculates, by the usual methods for statically indetermi- 
nate systems, the effect of the ribs in the case of unsymmetri- 
cal loading of the spars. The wing considered had nine like 
ribs. It was assumed that the middle rib was fixed, so that 

the calculation, simply concerns a wing structur e consisting of 

♦"Per Einfluss der S^ante"lraf~die Festigkeit der Holme." Tech- 
nische Berichte der Flagzeugmeisterei, Vol. Ill, (1917-1918), 
pp. 100-111. 
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two spars and four ribs* On the assumption that the same load 
was applied at every joint, Ball ens tedt calculates the course 
of the bending and torsional moments along the spars and shows 
the diminution of the maximum bending moment due to the union 
of the ribs # 

2. Dr. Thalau f s papers contain much valuable material.* 
He also discusses a rigidly supported cantilever wing, consist- 
ing of two spars and ribs. He defines this structure as a lat- 
tice girder. He chooses, as the expression of the reciprocal 
effect of the union of the spars, the diminution of the bending 
moment at the point of fixation: of the spars (fixation cross 
section) and shows that the end rib has the greatest effect and 
that the other ribs have only a negligible effect on the bend- 
ing moment at the point of fixation. Hence he investigates 
especially the shearing forces produced by the end rib in the 
case of like tapering spars*. 

3. Biezeno, Koch and Koning replace the rib union by a 
continuous union effect, as is done in the present paper.** 
They disregard the torsional rigidity of the ribs and consider 
them as inflexible. They give an approximation. method for the 

case when the rigiditi es of the two spars have the same ratio 

*"Zur Berechnung freitragender Flugzougflugel. « Z.F.1I., 1934, 
p.103; 1925, pp # 86-87 c "Uober die Verbundwirkung von Rirvpen 
im jf freitragenden, zweii holmigen und verspanungslosen Flugzeug- 
flugel, Z.F.M. , 1925, No. 20. 

**"Ueber die Berechnung von-, freitragenden Flugzeugf lugeln. " 
Zeitschrift fur Mech. und angew. Math. 1926, pp. 92-105. _ 
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in every cross section. Lastly, they investigate a wing with 
constant spar and rib cross sections on the assumption of uni- 
form loading and calculate the diminution in the bending moment 
due to the union. 

4. Reissner investigates the case of a wing situated 
above the fuselage and having two like spars, on the assumption 
that the stresses are constant in every cross sectiom of the 
spars (wings of constant strength).* Reissner obtains exceed- 
ingly simple formulas on the distribution, of the bending moment 
and on the reciprocal reduction of the load on the spars, for 
both symmetrical and unsymmetrical loading. For this case we 
can obtain, a simple expressions of the magnitude, which may be 
defined as the torsional rigidity of the whole wing. Reissner 
then employs this principle to determine the limits of static 
stability of wings exposed to the action of the air. 

5. Von Karman, in the winter semester of 1925/6 in his 
lectures on "Introduction to the Study of Airplanes," gave the 
general equations likewise on the assumption of an infinite 
number of ribs, but without the hypothesis of inflexible ribs 
and with consideration of the torsional rigidity of the ribs 
for a two-spar system and showed the influence lines for the 
loading of a spar by a single load applied to the other spar. 
We shall follow Von Karman' s lectures in the application of 

the differential equations. 

*"Heuere Probleme aus der Flugzeugstat ik, Verbundwirkung des 
Flugelkorpers." Z.F.M., 1926, pp. 181-185. 
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Bending and Torsion of a Cantilever Wing 

We investigated an unbraced cantilever wing structure, 
formed by two spars of uniform cross section and a number of 
equal ribs distributed at regular intervals along the spars. 
The structure was rigidly supported at one end. For the fol- 
lowing calculation, however, we assumed an infinite number of 
element ajy ribs regularly distributed along the spars and re- 
placing, in their effect, the ribs actually present. The ribs 
were assumed to be inflexible, but to have finite torsional 
rigidity with respect to their longitudinal axes. 

If both spars were stressed by separate or distributed 
loads proportional to the bending rigidities of the correspond- 
ing spars, the latter would be equally deflected. In this case 
the behavior of the spars was not affected by the ribs. 

Let p (x) and p 2 (x) denote the loads distributed 
along the spars; A x and A 2 , their corresponding bending 
rigidities; and x, the coordinates along the spars. If we 
then have the relation 

p x (x) p 2 (x) 

-a— '-IT- (1) 

the wing is simply bent. 3y a corresponding arrangement of 
the spars with respect to the wing chord, the above condition 
can be fulfilled in horizontal flight. This condition is not 
fulfilled in diving and gliding rfXTght, due to the travel of 
the center of pressure, and both spars are twisted about their 
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longitudinal axis due to their union through the ribs. In other 
words, they are stressed simultaneously with respect to both 
bending and torsion. From the assumption of absolute inflexi- 
bility of the ribs, it follows directly that the angles of 
torsion of both spars have the same value in every wing section. 
We call this angle the angle of torsion of the wing in the cor- 
responding cross sections. 

If the resultant of the external forces for all the wing 
sections lies between the two spars and divides the distance 
b (Fig. 1) into two parts b x and b 2 , so that b^bs = A 2 /A x , 
the wing then undergoes a simple bending. The point .0 is 
called the "elastic center" of the given wing section. The line 
connecting the elastic centers of all the wing sections is 
called the "elastic axis" of the wing* Obviously, the location 
of the elastic axis of a wing depends on the location of the 
spars with respect to the wing chord and on the bending rigidity 
of the spars. The elastic axis is a straight line in wings 
with constant or linearly diminishing bending rigidity. 

If we consider a wing with spars of constant rigidity and 
subjected to a load p distributed uniformly throughout their 
length, the wing is then subjected to a uniformly distributed 
load of 2 p and a moment of p (b x ~ "b s )« If we assume, on 
the other hand, that the spar cross sections, belonging to the 
same wing section, are subjected to two equal and opposite 
forces, the elastic axis will then suffer no deflection. This 
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is designated as the case of simple torsion. It is obvious 
that the general case can then be limited, by a corresponding 
shifting of the load, to these two individual cases. Since the 
calculation of the simple bending offers no difficulty, we can 
limit our calculations to the ce.se of simple torsion. 

The Torsional Rigidity of the Wing Section and 
of the Wing Itself 

While, for an elastic girder (generally of variable cross 
section) the ratio between the torsion per unit length and the 
torsional moment in each cross section; has a value dependent 
only on the cross section itself and characterizing the tor- 
sional rigidity of the section, this is not the case for a wing 
with spars of constant cross section. 

In fact, the torsional .angle & of a wing does not depend 
alone on the torsional moments in the two-spar sections, but 
also on the flexibility of the spars. The torsional moment is, 
in fact, equilibrated in every cross section in part by the tor 
sional moments of the spars and also by the moment of the shear 
ing forces in the spar sections. Obviously, these depend on 
the flexibility of the spars. Since the flexibility of the 
spars depends on the load distribution, it follows that the val 
ue & ! = d#/dx (x being the coordinate along the wing axis) 
in each wing section depends not only on the magnitude of the - 
torsional moment, but also on the load distribution. In other 
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words, the value £ f for each section of a cylindrical girder 
is proportional to the prevailing torsional moment, independ- 
ently of the value of the moment in the other sections of 
the girder. For a wing section., however, depends on the 
value and on the distribution of the torsional moment along the 
wing. 

For every wing section), we can therefore define by the 
equation 

Mt « ~ B£ f (2) 

a torsional rigidity B, which, however, depends both on: the 
characteristic quantities (bending and torsional rigidity) for 
the spars and ribs and also on the distribution of the torsion- 
al moment along the wing. Furthermore, by the expression:. 

Bf = M|i (3) 

we can introduce a torsional rigidity Bf for the whole wing. 
Thereby Mt is the mean torsional moment of the wing; l , its 
length; and & its distortion. 

We will undertake to determine the values B and Bf for 
a wing of constant spar and rib section in the following cases. 

1. Wing with a torsional moment increasing linearly from 
the tip toward the root. 

2. Wing with a torsional moment applied at the tip. 

In what follows it will be shown for both cases, to what 
extent the torsional moment in each wing section: is balanced 
by the torsion: of the spars and by the shearing forces in the 
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spar section* 

The Differential Equations of the Spar System on 
the Assumption of Inflexible Ribs 



We place the x-axis of the coordinate parallel to the 
spars with the origin at the wing tip, the direction from the 
wing tip toward the wing root being considered positive. The 
downward direction of the axes £ x and £ 2 is considered posi- 
tive (Fig. 2). The following symbols are used. 

B r> the torsional rigidity of a rib per unit length; 

n = l/a, the number of ribs per centimeter length of 
the spars; 

(3, the torsional rigidity of the ribs per centimeter 
length of the spars (p = n B r ); 

t x and £ 3J the bending deflections of the front and rear 
spars, considered positive downward (C = S x - £ 2 ); 

the torsional angle of any wing section with the ab- 
scissa x (i.e., the torsional angle of the spars 
with respect to the fixed root section, read posi- 
tively from the wing tip in the direction traveled 
by the hands of a clock); 

A I and A 3 , the respective torsional rigidities of the 
front and rear spars; 

b, the constant distance between the spars; 

\i x and |i 2 , the torsional moments per centimeter length, 
which the ribs transmit to the spars, these moments 
being considered positive in the direction of the 
motion of the hands of a clock as seen from the wing 
tip; 

t x and t 2 , the bending moments per centimeter length, 

which the ribs transmit "to the spars. The same num- 
bers with opposite signs represent the torsional mo- 
ments at the^ends of the ribs. r 1 and (applied 
to the spars) are considered positive, when they act 
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in the direction which conveys the positive £ axis 
in the positive direction of the x axis (Fig. 2); 

p x and p 2 , the external loading of the spars per centi- 
meter length; 

q x and , the loads per centimeter length, which the 

ribs exert on the spars due to the distortion: of the 
wing. 

a) The equilibr i um of the ribs .- The condition: of equilib,, 
rium for the forces (Fig* 3) reads q 1 + q a = 0 and for the 
moments 

q,- ^ (4) 
From these two equations we obtain 

% = _ (5) 



The condition of equilibrium for an elementary rib, with 
respect to rotation about its own axis, is 

t x + t 2 = 0, 

that is, 

T x = - T 2 - T (6) 

b) The equations for the bending of the spars »~ for any 
cross section; (abscissa x) of the front spar, the fundamental 
equation is 

d 3 L 
1 d x* 1 

in which VL X denotes the moment of the forces acting, between 
the wing tip and the given cross section, on the center of grav- 
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ity of said cross section* 

If we designate by Q >1 the shearing force in the given 
cross section: and consider it positive when directed upward 
(Fig. 2), we have 

d M, 



d x 



* = - Qj + T 



and consequently, 



Moreover, (Fig* 2) 



and hence 



=-Q, +t (7) 

1 & X 3 



A, fL-L = P x + 1, + £ (8) 
d x 4 1 1 dx 



For the rear spar, with the aid of equation (6), we like- 
wise obtain 3 

A 2 d \ = ~ Q s - t (7a) 
d x 3 

and 4 

2 d x 4 2 * dx 

o) The torsional equations for the spars e For the front 
spar the torsional equation* reads Bi £ f = - IL 19 in which M x 
denotes the torsional moment (about the center of gravity of 
the given cross section) of all the forces acting on the spar 
between this section and the wing tip. Taking into account 
the hypotheses concerning the positive direction of & , we 
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must consider M x positive when, as seen from the wing tip, 
it acta in the direction traveled by the hands of a clock. It 
is diq/dx = \x lf in which \x x is the torsional load of the 
front spar (and also the bending load of the ribs), considered 
positive when acting in the direction: of motion of the hands 
of a clock, as seen from the wing tip. We obtain 

Bl r - - *£. - - 

and, for the rear spar, B 2 a" = - |J S# The addition, of these 
two equations gives 

. (Bx + B 2 ) # M = — ( M*i + 

With the aid of this equation we can transform equations (4) 
and (5) into 

.q, = - 2i_122. (4a) 

and 

B, + B. 

% = V * (5a) 

d) Relation between the Flexure £ and the 
Torsional Angle £ of the Spars 

From Figure 3, it follows that 
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e) Relation "between the Flexure of the Spars and 
the Torsion of the Ribs 

Bearing in mind (Fig. 2) that the angle of torsion of a 
rib: of the length b is given by 

dx d x d x 

at every point of the abscissa x (cf. equation (9)), we ob- 
taim 

ids (10 ) 

b dx 

and, from this, 

dx b dx 



a (ID 



If we introduce the value of d^/dx into equation (8) 
and (8a), we obtain 



d 4 L " B d s £ 

On combining these two equations, we obtain 

Ilk *Lk 1 i+i 

d x* " d x 4 X 1 Pl ~ A 2 P a Q i ~ 

I q ,ifi *LVJL + iA 

A 2 2 b dx 2 VA X A 3 > 
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If we put i = -JL + JL and consider equations (4a), (5a), 
A A x A 2 

and (9), we obtain 

A t> s b U x aJ 

For simplification we put 

X ~ A b 8 

and 



K 2 = i + — ^ — = i + - 

B x + B 3 d B x + B 2 

and obtain 

a 1 r f 1 - k 2 k 2 *" + i (h, _ h>i = o, (12) 

as tlie- final result. 

Relation between the Shearing Forces and Torsional Moments 
in the Cross Sections of the Wing Spars 
Subjected to Simple Torsion 

The -torsional moment is held in equilibrium in every ling 
se<rtion -by *the shearing forces and torsional moments prevailing 
in the section. We will consider a wing subjected to simple 
torsion, being the torsional moment in any wing section* 

The condition of equilibrium for the portion of the wing includ- 
ed between the given section: and the wing tip (the forces and 
moments in this portion acting in the opposite direction to that 
indicated in Figure 4) is Q x = - Q 2 = Q for the forces and 



- -i a. b B r 
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-Qb - Mj - M 2 + Mt = 0 
for the moments. Since 

M x = - B x $' and M 3 = - B 2 

it follows that 

-Qb + (B 1 + B 2 ) a* + Mt = 0. (13) 
It now follows from equations (7) and (7a) that 

d x b dx d x 3 b dx 

and therefrom 

Q = - A x £ x ' " A 3 Ss' " £ f (14) 
2 " b ^ 

Ax C X ,M + A 2 £g = 0. (15) 

According to system of notation employed 

Al r' = a, (e 1 ,M - ) (is) 

A s C ,M = A 2 (e i M * - £ 2 '") (I?) 
From (15) and (16) we obtain by subtraction! 

-A, r = U x + A 2 K 2 ,M . 
From (15) and (l?) we obtain by addition 

A 2 r = U x + A 2 )£ 3 '" , 
from which it follows that 

* i 1 1 A 2 t-iit .frt A. c,tii 
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Multiplying by A x and A 2 we obtain: 

A, = A V" A s = - A £»", 

whence 

At" 1 = A i ^l'" " A s Sa* " 



Equation (14) then becomes 



Q = - AT" +| V. 
and, with the use of equation (9) 

Q = A b - £ * f 

and the condition of equilibrium (13) becomes 



*'" -w (Pb + Bi + BJ *' -ip 



or, with the use of the chosen abbreviation, 



£»" - K 3 \ 2 *' = . (18) 

From equation (13) we obtain 

This equation shows the relation between the shearing forces 
and torsional moments of the spajs in every wing section. 
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Case I. Spars with Uniformly Distributed Load 
(p x and p s being constants) 

The fundamental equation (12) holds good for any distribu- 
tion of the loads p x and along the spars, hence also for 
the case when a single force acts on the wing tip. We will con- 
sider the case when the spar loading is uniformly distributed, 

1 /P P \ 

i.e., when p, and p P are constants. The member — ( — L 2-) 

12 bUi As/ 

of the differential equation is then a constant. The general 
equation for this case is 

-KXx 

£ = C 1 e + C 2 e + C 3 x -f 

C 4 + 3_ (*l -2^x . 

2b x 2 \ 2 V Ai A2 y 

in which the four integration constants C iy C 2 , C 3 , and C 4 
are to be determined from the following limiting conditions. 
For x = I we have, as a result of the rigid fixation, S- = 0 
and = 0, instead of x = 0, that is, at the wing tip we 
must have £ n = 0, due to the vanishing bending moment £ tf . 
If it is noted that the torsional moment is zero at the wing tip, 
it follows for equation (18) that 

1 f - K 2 X s £ r = 0 

x=o x=o 

is the fourth limiting condition. 
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For these limiting conditions the constants "become 
c = - 1 fl_ 

1 2b K 3 \ 3 \t ' " 2 ^ 



2b K X 3 VAj. A 2 / K z X cosh K X I 

0 1 / Pl M Hi - e K ^ 

2 2b K 2 X 2 ^A, A 2 ' K 2 X 2 cosh K X I 

C 3 = 0 

n = (2± - ^ — iK 2 X 2 I 3 cosh KXl- 

2b K X 2 VAj. A 3 / K 2 X 2 cosh K X I 1 



-2 K X I sinh K X I - 2 \, 



and we have 



Pi. _ Ps. 

£ = _ Ai Ai J. 2 cosh KX (I - x) - 2 K X I (sinh KXl- 

2b K* X 4 cosh K X I I 

- sinh K X x) + K 2 X 2 (l 3 - x 3 ) cosh K X l - a}' 

When = & then * = 0. 
A x A a * 

Hence we can always resolve the case when p x and p 2 have 
any desired values into two cases, that of simple bending, in 
which the load is proportional to the bending rigidity, and that 
of simple torsiom, in which the loads of the two spars are alike, 
and opposite. 

For the latter case, we obtain the formula 

From here on we will confine our attention to this case. The 
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load can "be so expressed that the wing will be subjected to a 
linearly variable moment of the magnitude Mt = p b x. 

a) The torsional angle is expressed by the equation 

£ = — ■ / 2 cosh K a. ( l - x) - 

2b A K X cosh K X I 1 

-2 K X I (sinh K X l - sinh K X x) + 

+ K 2 \ 2 (I 2 - x 2 ) cosh K X I - 2 } . 

Since 

b 2 A X 2 = B x + B,, 

we can write: 

* = P_L_i! 1 /a cosh K X (l - x) 

2 X 2 (B,. + Bg) K s X 2 I s cosh K X I 1 

-3KH (sinh K X I - sinh K X x) + 

+ K 2 X s ( I s - x 2 ) cosh K X I - 2) (20) 

If we put Bi = K 3 (B x + B 2 ), in which Bi may be regarded 
as the sum of the torsional rigidities of the spars and of the 
ribs situated within the distance b, then 

* = p b 1 — -2 h cosh K X (I - x) - 

2 Bi K X 2 I cosh n M 

-2 K X I (sinh K X I - sinh K X x) + 
+ K 2 X 2 (I 2 - x 2 ) cosh K X I - 2 }. 

The term p b l 2 /2 Bi represents the angle of torsion 
■&io of a girder of the torsional rigidity Bi. We then obtain 
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■J— = — ; rs-^- {z cosh K X (I - x) - 

2 K X I (sinh KXl- sinh K X x) + 

+ K s X 3 ( I s - x 2 ) cosh K X I - 2 } (21) 

*/*io is plotted against x/l for different values of 
K X I in Figure 5 (Table I). 



Values of: 
1 



10 



(KXl) 2 cosh K X I I 



TABLE rfelg. 5) 



«2 cosh KX(l-x) -2 KXl (sinh KXl 



sinh KXx) + K 2 X s (I 5 -x s ) cosh K X I - 2 } 



x/l 


KXl=l 


KXl=2 


KXl =3 


KXl=4 


0.0 


0.181 


0.403 


0.530 


0.630 


0.2 


0.134 


0.314 


0.440 


0.537 


0.4 


0.888 


0.217 


0.325 


0.404 


0.6 


0.046 


0.124 


0.191 


0.247 


0.8 


0.014 


0.038 


0.057 


0.087 


1.0 


0.000 


0.000 


0.000 


0.000 



From equation (21) we obtain for x = 0 



10 



(KXl) 



2 tanh KXl 
K X I 



(KXl) 3 cosh KXl 



As KXl approaches » , this ratio approaches 1, which means 
that the wing undergoes a torsion equal to that of a girder of 
the sajne length with the torsional rigidity B±. 

^o/^ix) is plotted against KXl in Figure 6 (Table II). 
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TABLE II (Fig. 6) 
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Values of: 

10 



= 1 + 



(KXl) ; 



3 tanh KM 
KM 



(KM) cosh KXl 



KM 


*o/*io 


KM 


V*lo 


KXl 


V»to 


0.0 


0.000 


0.8 


0.132 


8.0 


0.779 


0.1 


0.002 


0.9 


0.153 


10.0 


0.820 


0.2 


0.010 


1.0 


0.181 


12.0 


0.847 


0.3 


0,023 


2„0 


0.407 


16,0 


0.882 


0.4 


0.040 


3„0 


0.539 


20.0 


0.905 


0.5 


0.0-32 


4-0 


O.o30 


25.0 


0.923 


0.6 


0.080 


5.0 


0*830 


30.0 


0.936 


0-7 


0.105 


6.0 


0.7*2 








Torsional 


rigidity of the wing.- 


According to 


the pre- 



vious definition, the torsional rigidity of the wing is 

Bf = ^ £ 0 • For the win S "tip (x = 0) we obtain from equa- 
tion (20)* 



tanh K X I 



0 K 2 (B 1 + B a ) '< K 3 X s I s K X I 

+ * 



Consequently 
Bf = 



K 2 X s I 2 cosh K X l J 
K 2 (Bj + B 2 ) 



V 



1 + 



2 tanh K X l 



(K X if K X I (K X if cosh K X I 

As K X I approaches 0, Bf approaches «> and as K X I ap- 
proaches co, Bf approaches B±. For ribs without torsional 
rigidity, the torsional rigidity of the wing therefore coincides 
with the sum of the torsional rigidities of the spars. 
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However, if the ribs had absolute torsional rigidity, we 
would have Bf = «» for an infinitely long wing, i.e., the wing 
would not be distorted. Bf/Bi is plotted against K X I in 
Figure 7 (Table III). 

TABLE III (Fig. 7) 

Values of: 



3 f _ 1 

Bi ~ 1 + 2 _ 2 tanh KXl ~ 2 

(KXl) 3 KXl " (KXl) s cosh KXl 



KM 


Bf 

BT 


KXl 


Bf 

BT 


KXl 


Bf 

BT 


0.0 


oo 


1.0 


5.46 


4.5 


1.56 


0,1 


400,00 


1.2 


4.52 


5.0 


1.46 


0.2 


100.00 


1.4 


3.34 


5.5 


1.22 


0.3 


47.60 


1.6 


3.08 


6.0 


1.38 


0.4 


25.00 


1.8 


2.70 


8.0 


1.28 


0.5 


19.20 


2.0 


2.46 


10.0 


1.22 


0.6 


12.50 


2.5 


2.08 


12.0 


1.18 


0.7 


9.54 


3.0 


1.86 


16.0 


1.12 


0.8 


7.68 


3.5 


1.70 


20.0 


1.10 


0.9 


6.54 


4.0 


1.60 


30.0 


1.06 



c) Torsional rigidity of the wing sections . ~ This is 
B = - ^ = p (see equation (2) ). is obtained from 

equation (21) by differentiation according to x. 

a* = - — — ^- . .. „ . . Jz X (I cosh K X x - x cosh 2 X I) 

Bi K X cosh K X I I 

- sinh K \ (I - x) } (22) 

Herewith we obtain 

B _ K X x cosh K X I 

Bi K X ( I cosh K X x - x cosh KXl)- sinh K X (I - x)' 
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In Figure 8 (Table IV) B/B\ -is plotted against K/l with 
K X I as parameter. B/Bj_ = » when x = I, due to the assumed 
rigidity of fixation of the wing root. 



TABLE IV (Fig. 8) 



Values of: 



B _ K X x cosh K X I 

Bi KXTl cosh K X x - x cosh K X I) - sinh K X (I - x) 



X 

I 


XX 1=1 


KX 1=2 


KXl=3 


0.0 


0.000 


0.000 


0.000 


0.2 


1.748 


0.876 


0.762 


0.4 


2.570 


1.632 


1.240 


0.6 


5.6*20 


2.536 


1.752 


0.8 


12.488 


4.710 


2.980 


0.9 


42.200 


9.000 


4.320 


1.0 


00 


oo 


CO 



d) Course of the torsional moments along the spars .- The 
formulas B x £ f and B 2 give the torsional moments transmit- 
ted "by a single spar section. Table V gives us the values of 

, that is, the angle of torsion per unit length, with respect 
to the mean angle of torsion, ify' = - p b l/B±. We have 

*' = 2E _ K X 1 cosh K X x - sinh K X (I - x) 
~$±~ I K X I cosh K X I 

In Figure 9 ty^i' is plotted against x/l for different 
values of K X I (Table V). 
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Values of: 



x K X I cosh K X x ~ sinh K X ( I - x) 





I 


K X I cosh K X I 




X 
I 


KXl=l 


KXl=2 


KXl=3 


KXl=4 


0.0 


0.113 


0.216 


0.232 


0*250 


0.2 


0„114 


0.228 


0.281 


0.280 


0.4 


G.113 


0.245 


0.317 


0.355 


0.6 


0-097 


0.236 


0.331 


0.418 


0.8 


Co 06 4 


0.170 


0.270 


0„358 


1.0 


0,000 


0.000 


0.000 


0.000 



The table and curves can be used for plotting the torsional 
stress of the spars. The maximum torsional stress does not 
come at the ends of the spars. 

e) Course of the shearing forces in the spar sections .*- 
In this case equation (19) reads 



Bt "J - Bp _ , t 

— 5- £ f + p x 



or 



Q = — u 

_g_ . x ; + b. >t- 

p x b 



Introducing £ ! from equation. (22) we obtain 

Q s= i + — J Jl cosh K X x sinh K X (l - x) _ i \ (23) 
p^x K s lx cosh K X I ~ K X x cosh K X I J 

For every spar section, the ratio Q/px determines the portion 
of the external shearing forces held in equilibrium by the tan- 
gential stresses in the cross section. For x = I, Q/p^ = 1> 
that is, the whole external load, which is applied to the spar, 
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is balanced "by the tangential stresses in the root section* 
This indicates that 

I 

Qx=o + / q d x = o, 

0 result 
i.e., the resultant of the forces acting on the spars, as the/ 

of their being connected by the ribs, is zero.* In other 
words, in a rigidly mounted wing the spars cannot mutually re- 
lieve one another, Q/px is plotted against x/l in Figure 10, 
for different values" of K X I and K (Table VI). 

TABLE VI (Fig. 10) 

Values of: 



ft -1 + 1/1 cosh K X x _ sinh K X (I - x) _ ^ 
p~ x ~ K 2 ^x cosh K X I K X x cosh K X i S 



X 

I 


k n = i 


n t = 2 


K=l 


K=1.2 


K=1.4 


K=l 


K=1.2 


K=1.4 


0.0 




— OO 


— 00 


— OO 


— 00 


— 00 


0.1 


-0.15 


0.303 


0.416 


-6.83 


-4.44 


-3.00 


0.3 


0.45 


0.618 


0.730 


-3.41 


-1.37 


-0.74 


0.3 


0.62 


0.739 


0.808 


0.31 


0.44 


0.60 


0.4 


0.73 


0.806 


0.857 


0.39 


0.57 


0.69 


.0.5 


0.78 


0.835 


0.890 


0.51 


0.65 


0.75 


0.6 


0.83 


0.883 


0.914 


0.50 


0.73 


0.79 


0.8 


0.93 


0.946 


0.960 


0.79 


0.85 


0.89 


1.0 


1.00 


1.000 


1.000 


1.00 


1.00 


1.00 



f ) Course of the load reaction q along the spars *- For 
every spar section at the distance x from the wing tip, we have 
the equilibrium condition 

x 

p x ~ Q = Q x =o + / q d x 

♦First discovered by Biezeno, Koc^i and Koning ("Ueber die Berech- 
nung von freitragenden Flugzeugflugeln. " Zeitsch. f. angew # 
Math. u. Meoh. , 1926, pp. 97-105.) 
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According to equation (22) 

Q 1(1 cosh K X x sinh K X ( I - x) 



1 - 

p x 



/ L cosh K X x sinh K X ( I ~ xj ^ "\ 
K 3 \x cosh K X I ~ K X x cosh K X I f 



and hence 



r* q / cosh K X x sinh K X ( I - x) ^ 

Qx=0 V Q * " " F I 1 cosh K X I ~ K X cosh K X I *i ; 

By derivation according to x 

q = ~IL-/x\ i sinh K X x + cosh K X ( I - x) _ 

K 2 I ' cosh K X I ~ cosh K X I J 

For x = 0 (wing tip), q = 0; For x = I (wing root), 

q = ^ P /k X I tanh K X I + — ll. 

K 2 I cosh K X I J 

Figure 11 shows the curve of 

K 2 Q = K X I sinh K X x + cosh K X ( I - x) 



cosh K X i 

plotted against x/l for different values of KM. See Table 
VII for the computed values. As shown by Figure 11, q indi- 
cates an increase in the load p toward the wing tip and a de- 
crease toward the wing root. 
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Values of: 

v s q _ K X I sinh K X x + cos] 
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p 




cosh K X I 




T 




i 1 ■" 
KXl=2 


KXl=3 


0.0 


0.000 


0.000 


0.000 


0.2 


0.010 


0.090 


0.259 


0.4 


-0.034 


0.047 


0.242 


0.6 


-0.112 


-0.157 


-0.055 


0.8 


-0.237 


-0.547 


-0.745 


1.0 


-0.410 


-1.191 


-2.099 



g) Dependence of force Q n on K X I .- For x = 0, equa- 
tion (19) gives 

0 b 0 

Introducing x = 0 into equation (22) gives 

p b I 



K 2 (B x + ) cosh KUl 
whence we obtain 



j sinh K X I 



K X I 



p 



2 - 1 ( 1 _ 

I K s \ cosh K ^ I 



tanh K X I 



K X I / 



As K X I approaches » , ft 0 /p I approaches 0. Q 0 /p^ is 
plotted against K X I in Figure 12. The curve reaches its max- 
imum value at about K X I = 3. 
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Values of: 



pi 



x 2 ( 



TABLE VIII (Fig. 12) 

1 tanh K X I 



cosh K X I 



KX I 



) 



KXl 


K=l 


K=1.2 


K=1.4 


0 


OoOOO 


0.000 


0.000 


1 


-0< 113 


-0.077 


-0.057 


2 


-0.216 


-0.150 


-0.110 


3 


-0,233 


-0.165 


-0.118 


4 


-C.213 


-0.148 


-0.108 


5 


-0.186 


-0.129 


-0.095 


6 


-0.162 


-0.112 


-0.082 


8 


-0.124 


-0.086 


-0.063 


10 


-0.100 


-0.069 


-0.051 



Case II. T/Ting Subjected to a Constant Torsional Moment 
Applied to Its Tip 

We shall assume that the wing is subjected to the action 
of a torsional moment, exerted, e.g., by two equal and opposite 
forces acting on the spars. These produce the same torsional 
moment Mfc = P b, in every wing section. Since p x = pi, =0, 
the fundamental equatiom (12) now becomes 

The general integral of this equation is 

KXx -K\x 



£ = 0 1 e + C 2 e 



+ C 3 x + C 4 , 



where C x , Cg , C3 and C 4 are the four integration constants, 
The limiting conditions for this case are 
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^x= 1=0 
*'x=l = 0 



4"x=o = 0 



whence we obtain 

M 4 



o s = - 



2 A b 2 K 3 X 3 cosh K X I 
M t 



2 A b 2 K 3 X 3 cosh K X I 



A b 2 K 2 X 2 



_ M-fc I / tanh K X I ^ 

4 " K 2 X 2 A b 2 V K X I 1 ) 



and 



M+ 

* = -= — 5 5— sinh K X x - sinh K X I + 

K 3 X 3 A b 2 cosh K X I 

+ K \ (l - x) cosh K X I (24) 

a) Course of the torsional angle along the win g.- If we 

introduce X" = ^ 1 + f* 2 as in equation (12), we can then write 

A b~ 

equation (24) in the form 

& = — : 5 J sinh K X x - 

Kf (B x + B 2 ) K ' X I cosh K X I 1 

- sinh K X I + K X I (I - x) cosh K X I j- (24a) 
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According to St, Venant ! s theory 



*io = 



M t I 



K 2 (B x + B a ) 

indicates the maximum angle of torsion of a girder of length I and 
torsional rigidity Bi = K 3 (B 1 + B a ) • Then 

* _ l 



sinh K X x - 



'io KXl cosh K X 

- sinh KXl + K\(l-x) cosh KXl} (25) 

In Figure 13, */*io is plotted against x/l for various con- 
stant values of KXl. The dash-dot line here represents the 
function &i/$i 0 y in which $± is the angle of torsion along 
the girder, 

TABLE IX (Fig- 13) 

Values of: 

£ _ „ x tanh K X I 



'10 



= i - * 
i 



K X I 



sinh K X x 
K X I cosh K X I 



X 

I 


KXl=l 


KXl =2 


KXl =3 


0.0 


0.240 


0.518 


0.672 


0.2 


0.170 


0.373 


0.490 


. 0.4 


0.105 


0,236 


0.318 


0.5 


0.050 


0.118 


0.167 


0.8 


0.014 


0.032 


0.049 


1.0 


0.000 


0.000 


0.000 



From equation (25) we obtain for x = 0 



ip_ _ i tanh K X I 



no 



K X I 
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^o/^io is plotted against K X I in Figure 14 (Table X). This 
ratio approaches 1, as K X I approaches <*> . i n this case the 
wing behaves like an ordinary girder with a torsional rigidity 
equal to the sum of the torsional rigidities of the spars and 
ribs. 

TABLE X (Fig. 14) 

Values of: 



= 1 ~ 



tanh K X I 



KXl 



KXl 


"^io" 


KXl 


Ho 


KXl 




0.0 


0.000 


1.0 


0.240 


4.5 


0.779 


0,1 


0.003 


1.2 


0.306 


5.0 


0.800 


0.2 


0.014 


1.4 


0.367 


5.5 


0.818 


0.3 


0.024 


1.6 


0.426 


6.0 


0.835 


0.4 


0.050 


1.8 


0.473 






0.5 


0.075 


2.0 


0.515 






0.6 


0.105 


2.5 


0.604 






0.7 


0.136 


3.0 


0.679 






0.8 


0.170 


3.5 


0.714 






0.9 


0.205 


4.0 


0.749 







M+ 1 

b) Torsional rigidity of the win g.- This is Bf = . 
For x'= 6, equation (24a) gives 



/, Mt I- r tanh K X I -\ 

~ + 3=) "I win r 



o - 



Hence 



K 2 (B x + B 2 ) 



Bf 



K X I 



B; 



n tanh KXl 
1 K X I 



Bf/Bi is plotted against KXl in Figure 15 (Table XI). 
Bf approaches », as KXl approaches 0. 
Bf approaches Bi, as KXl approaches ». 
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TABLE XI (Fig. 15) 

Values of: 



Bf = 1 

Bi tanh K X I 

KXl 



KXl 


Bf 
Bi 


KXl 


Bf 
Bi 


KXl 


Bf 
Bi 


0.0 


00 


1.0 


4.18 


4.5 


1.28 


0.1 


303.00 


1.2 


3.27 


5.0 


1.25 


0.2 


7.1.50 


1.4 


2.72 


5.5 


1.22 


0.3 


35.70 


1.6 


2. 35 


5.0 


1.20 


0.4 


20.00 


1.8 


2.11 






0.5 


13.20 


2.0 


1.94 






0.6 


9.52 


♦2.5 


1.65 






0.7 


7.36 


3.0 


1.49 






0.8 


5.90 


3.5 


1.40 






0.9 


4.88 


4.0 


1.33 







o) Torsional rigidity of the wing section .- This is 
B = - lit/*' (cf. equation (2) ). When differentiated with re- 
spect to x (equation (24a) ), 

£»= !!* cosh K X x - cosh K X l\ (26) 

K 2 (B x + B 2 ) cosh KXl J 

Hence 

B = K 2 (B, + BgJ ^ 
1 - cosh K X x ' 
cosh KXl 



K 2 (B 1 + B 3 ) ~ Bi ' 1 _ cosh K X x * 

cosh KXl 

At the wing root, due to the assumed rigid fixation, we have 
B/Bi =«. B/Bi is plotted against xA for different values 
of KXl in Figure 16 (Table XII). 
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TABLE XII (Fig. 16) 

Values of: 

B 1 



Bi ^ cosh K X x 
cosh K X I 



X 

I 


K\l=l 


K\l=2 


K\l=3 


0.0 


2.840 


1.360 


1.115 


0.2 


2.950 


1.400 


1.142 


0.4 


3.350 


1.555 


1.230 


0.6 


4.320 


1.920 


1.445 


0.8 


7.500 


3.18 


2.230 


0.9 


14.300 


5.75 


3.780 


1.0 


OO 


oo 


00 



a) Course of the shearing forces in the spar sections alon g 
the win g»~ In this case equation (19) reads 

Q = h^JjL **' + p. 

Using equation (36) for we obtaini 

Q = ? , J C osh K X x - cosh K X I X + P 

K cosh K \ M / 

and 

Q _ 2, . , 1 (l - cosh K X x^ 

p - ~ K 3 V cosh K X 1/ ~ 

For x = I (wing root) we- obtain?. Q/p = 1. For x = 0 (wing 
tip) we obtain 

So = i _ L- fi - i > 

P ~ K 2 V cosh K X I / * 

For very long or very flexible wings this equation can be writ- 
ten 

Q _ , 1 
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If K = 1, that is, if the torsional rigidity of the ribs is 0, 
we obtain Q/P = 0 for the wing tip, i.e., Q = 0. This means 
that the torsional moment of the wing near its tip is balanced 
only by the torsion of the spars. 

In the case when the ribs have an infinite torsional rigid- 
ity, Q = P for every spar section. The torsional moment is 
then balanced onl^ by the shearing forces acting in the spax 
sections. The shearing force Q 0 at the wing tip for x = 0 
(equation (27) ) is 

0 r V cosh K X I S 

Q 0 approaches P fl - as K X I approaches « . 

Qo approaches P, as K approaches » . 

Q 0 /P is plotted against K X I in Figure 17 (Tahle XIII ). 

TABLE XIII (Fig. 17) 

Values of: 



3fi = i _1 f l ~ 1 . > 
p K s V cosh K X Ij 



KXl 


K=l 


K=1.2 


K=1.4 


0 


1.000 


1.000 


1.000 


1 


0.648 


0.756 


0.821 


2 


0.2G3 


0.490 


0.626 


3 


0.099 


0.375 


0.540 


4 


0.036 


0.330 


0.510 


5 


0.013 


0.314 


0.496 


6 


0.009 


0.310 


0.493 


8 


~ 0.000 


0.307 


0.491 


10 


~ 0.000 


0.306 


0.490 
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From these considerations it follows that an increase in 
the torsional rigidity of the ribs entails an increase in the 
torsional rigidity of the wing and consequently a diminution of 
the torsional moment accompanied "by a simultaneous increase of 
the shearing forces in the spars. Q/P is plotted against x/l 
at different values of K \ I and K in Figure 18 (Table XIV). 

TABLE XIV (Fig. 18) 

Values of: 

9t = i _ _L h _ cosh K A, x \ 
P k 2 V cosh K \ I j 



X 

I 


K 


X I = 


1 


K 


X I = , 


3 


K 


\ l = 


3 


K=l 


K=1.2 


K=1.4 


K=l 


K=1.2 


K-1.4 


K=l 


K=1.2 


K=1.4 


0.0 


0.648 


0.756 


0.821 


0.2S6 


0.490 


0.626 


0.099 


0.375 


0.520 


0.2 


0.661 


0.765 


0.827 


0.283 


0.505 


0.637 


0.117 


0.387 


0.550 


0.4 


0.700 


0.792 


0.847 


0.358 


0.555 


0.671 


0.180 


0.330 


0.581 


0.6 


0.768 


0.839 


0.881 


0.487 


0.648 


0.735 


0.307 


0.519 


0.646 


0.8 


0.868 


0.908 


0.932 


0.686 


0.782 


0.840 


0.550 


0.687 


0.770 


1.0 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


l.oob 


1.000 



e) Cou rse of the supporting reaction q. along the wi ng • - 
It was found that equation (27) which, after transformation, reads 



cosh K \ x 



cosh K X I/* 

holds good for every spar section, Q is that portion of P 

which is balanced in every spar section by the tangential stresses 

The other part of p, namely, P -i ( 1 - cos £ £-v-r ^ is bal- 

lf \ cosh K \ V ; 

anced by the forces transmitted by the ribs to the portion of 
the spar between the wing tip and the given cross section. 
Therefore we may write 

Qo + A d z = - P \ (l ~ 
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by differentiating which we obtain 

M-t X sinh K X x 



K 2 cosh n I ' 
q = 0 for x = 0 (wing tip), 
q = tanh K X I for x = I (wing root) 

q approaches 0 as K approaches », that is, for 
ribs with absolute torsional rigidity, 



ql£ is 



X M t 

plotted against x/l for different values of K X I 
in Figure 19 (Table XV). 

TABLE XV (Fig. 19) 



Values of: 



b K 2 



sinh K X x 



XM-t cosh K X I 



X 

I 


KXl=l 


KXl=2 


XX 1=3 


0.0 


0.000 


0.000 


0.000 


0.2 


0.130 


0.109 


0.063 


0.4 


0.266 


0.236 


0.150 


0.6 


0.412 


0.402 


0.292 


0.8 


0.576 


0.632 


0.542 


1.0 


0.762 


0.964 


0.996 



f) Course of the torsional moment along the spars .- Equa- 
tion (36) was transformed into 



K a (Bi + B a ) _ cosh K X x _ -j_ 
Mt cosh XXI 



M t 

According to St. Venant's theory, the ratio 

M t = Mt 

X 2 (B x + B s ) Bi 



(28) 



N.A.C.A. Technical Memorandum Ho. 520 37 

represents the constant increment of the torsional angle of . a 
cylindrical girder of torsional rigidity B±* We put 



with which equation (28) can "be written 

_ -i cosh K X x 
W cosh K X I * 

TABLE XVI (Fig. 20) 

Values of: 



£ f = -i cosh K X x 
"V" ~ cosh K X I 



I 


K\l=l 


K\l=2 


K\l=3 


KXl=4 


0.0 


0.352 


0.735 


0.896 


0.964 


0.2 


0.339 


0.714 


0.876 


0.950 


0.4 


0.299 


0.644 


0.820 


0.906 


0.6 


0.232 


0.520 


0.692 


0.796 


0.8 


0.133 


0.315 


0.450 


0.550 


1.0 


0.000 


0.000 


0.000 


0.000 



This ratio is evidently proportional to the torsional mo- 
ment in every section, -i* 1 /^ 1 is plotted against x/l for 
different values of U I in Figure 20 (Table XVI). 
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S u m m a r y 



The present paper is confined to the case of a cantilever 
wing with constant spar and rib sections. Although this assump- 
tion is seldom borne out in practice, nevertheless the investi- 
gation of this simple case enables us to gain a clearer idea of 
how the distortion of a two-spar wing is effected. 

The quantity- 



first enters into the computation as the essential parameter. 
The tables bhA figures show that, for small values of this pro- 
portionality factor, the distortion depends chiefly on the bend 
ing rigidity of the spars. In other words, the distortion is 
the result of the difference in the deflection of the two spars 
At large values of XI, however, the wing behaves essen- 
tially like a twisted girder whose torsional rigidity equals 
the sum of the torsional rigidities of the spars (strengthened 
to a certain extent by the torsional .rigidity of the ribs)* 
The second parameter 



as previously remarked, expresses the effect of the torsional 
rigidity of the ribs. The tables furnish accurate information 
regarding this effect. For large values of X I, this effect 
can be so expressed that it is only necessary to add to the tor 
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sional rigidity of the spars the torsional rigidity of the ribs 
per length unit b (distance between the spars). 

All the calculations were made on the assumption that the 
system of ribs can be replaced by a uniformly distributed elas- 
tic union. It is a further task for theory and experimentation 
to determine to what extent this assumption will give sufficient- 
ly accurate results for practical purposes. 

Translation by Dwight M. Miner, 
National Advisory Committee 
for Aeronautics. 
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Figs. 8 & 9 
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Fig. 12 
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Fig. 14 
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Figs. 15 & 16 
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Figs. 17 & 18 
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